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We discuss memory effects in the conductance of hopping insulators due to slow rearrangements of 
many-electron clusters leading to formation of polarons close to the electron hopping sites. An abrupt 
change in the gate voltage and corresponding shift of the chemical potential change populations of 
the hopping sites, which then slowly relax due to rearrangements of the clusters. As a result, 
the density of hopping states becomes time dependent on a scale relevant to rearrangement of the 
structural defects leading to the excess time dependent conductivity. 
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Introduction - Memory effects in low-temperature 
transport properties of hopping insulators have been re- 
ported in several systems 0, 0, 0] ■ After excitation 
from equilibrium by, e. g., a sudden change of a gate 
voltage, V g , the conductance of the system increases in- 
dependent of the sign of the change. This excess conduc- 
tance, Act, may persist for long times after the excitation 
forming the so-called memory cusp, see [H for a review. 

Several concepts were used to explain the memory 
cusps in the dependence of the conductance, G, versus 
the gate voltage, V g . So-called intrinsic mechanism, is 
based on the assumption that the memory effects are 
due to complex dynamics in the strongly correlated sys- 
tem of interacting electrons @, 0, Q ■ It is a natural as- 
sumption since hopping insulators lack strong metallic 
screening and the long-range Coulomb interaction can be 
decisive. In [§] the connection between the glass-like be- 
havior and Coulomb gap was argued. Another scenario, 
so-called the extrinsic, assumes that the observed con- 
ductance relaxations are due to the influence of slowly 
relaxing atomic configurations acting on the conducting 
channels. It was first advocated in [hJ to explain the 
occurrence of a G(V g ) cusp in granular gold films. 

To the best of our knowledge, the mechanism behind 
the memory effects in hopping insulators is far from be- 
ing fully understood. In particular, we are not aware 
of fully convincing explanation of the experimental re- 
sults [j| on "double-dip" structure of the G{V g ) depen- 



dences, and their relaxations. Recent experiments [11 1 
aimed at studies of the influence of the sample lateral 
dimensions on the glassy properties show that there are 
reproducible conductance fluctuations having apparently 
different time scale comparing with the memory cusps. 

Recently we suggested a simple extrinsic model allow- 
ing for the "two-dip" behavior of the conductance of a 
structurally disordered hopping insulator [13]. Accord- 
ing to this model, the memory is supported by two-state 
dynamic structural defects present in any medium with 



sufficient amount of structural disorder [l3j. The two- 
level defects get polarized by the electrons and, in turn, 
form a polaron gap at the hopping sites decreasing hop- 
ping conductance. The slow dynamics of conductance is 
then due to slow rearrangement of polaron clouds around 
the hopping sites. An important feature of this extrin- 
sic model qualitatively explaining logarithmic relaxation 
and memory effects [3|, |4| is presence of a set of fluc- 
tuators possessing dipole moments, uniform density of 
states and logarithmically uniform spectrum of relax- 
ation times. Such fluctuators can also have intrinsic, 
electronic nature. Recent experiments [Tlj can help 
to determine which particular mechanism - "extrinsic" 
or "intrinsic" - is responsible for the observed behavior 
using the temperature dependence of the "typical" relax- 
ation time defined by the memory deeps equilibration. 
We will show that the "chessboard" electronic fluctua- 
tors suggested in [IH to interpret 1/f noise in hopping 
conductivity also lead to the non-equilibrium conductiv- 
ity behavior similar to that observed experimentally. As 
we will see below, this time is sensitive to the minor devi- 
ations of the relaxation time distribution from its cx 1/r 
shape. The model [l5[ turns out to be able to explain 
the observed increase of the typical relaxation time with 
temperature, as well as some other observed features. 

Main concepts - To calculate the time-dependent 
conductance of the system one needs to know the time- 
dependent density of hopping states (DOS) at the Fermi 
level. The DOS is modified by rearrangement of the pop- 
ulations of the sites neighboring the sites belonging to the 
hopping cluster. This rearrangement driven by hopping 
electrons lead to decrease, U, of the energy of a hop- 
ping site. The polaron shifts U are different for different 
sites and can be characterized by a distribution function, 
F{U, t) . This distribution slowly depends on time due to 
rearrangement of the polaron clouds. 

The analysis below is similar to our previous work [l2| . 
A site can be brought from the ground to an excited 
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state by placing or removing an electron, that can take 
place only if the excitation energy, s, exceeds the polaron 
gap. Thus the formation of the polaron gap excludes all 
the states with U > e. Here e is the typical excita- 
tion its equilibrium value being the width of the hopping 
band, Eh- For an equilibrium state created by sudden 
change SV g » Eh of the gate voltage e w SV g . Thus, 
e ~ max{e^, SV g }. The effect of the polaron cloud on the 
conductivity can be estimated as 

SG(t) 



G 



T{U, t) dU . 



(1) 



This result is straightforward - the sites inside the po- 
laron cloud cannot be occupied, and the density of states 
is less than the universal one. Thus the relative decrease 
in the density of states due to polarons can be estimated 
as a relative volume occupied by the polaron cloud. 

To find !F(U,t) one need to specify the slow relax- 
ing aggregates producing the polaron shift U. Follow- 
ing previous studies of 1/f noise in hopping conductiv- 
ity [15| we assume that slow dynamics is due to chess- 
board electro-neutral clusters having 2N sites placed at 
a distance s» r between each other. The relaxation rate, 
t of a given aggregate depends both on the number of 
sites, 2N, in the aggregate and on the typical distance, 
r, between the sites, and energy difference E between its 
lowest energy levels. We are interested in the domain 
where E < T since in the opposite case the cluster re- 
sides in its ground state. Let us define the distribution 
P(N, r, E) such that P(N, r, E) dN dr dE is the number 
of the clusters per unit volume having the parameters 
within the region (N + dN, r + dr, E + dE) . To estimate 
this distribution, we will take into account that the typi- 
cal aggregate volume is Nr 3 while the energy bandwidth 
for small energies is equal to ^/~Ne 2 /nr. Thus, the to- 
tal density of states, W(N,r,E), of all aggregates with 
the intersite distances r' larger than some r is given as 
W{N,p,E) ~ X N /T Q N 3 / 2 a 3 p 2 . Here A is a probability 
to add additional pair of sites to the aggregate; p = r/a 



where a is the localization length, while To 
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In this way we get the expression for the partial den- 
sity as P{N,p) ~ \ N /T N 3 ' 2 a 3 p 2 . Since we are inter- 
ested in the case E < T -C T the partial density is 
^-independent. Let us first consider the Efros-Shklovskii 
(ES) regime [3] of variable range hopping (VRH). In this 
case the relaxation rate for cluster rearr ang ement can be 
expressed by the interpolation formula [151 ] 



- 1 (N,p) = Vo (e~ NV ^'P + e- N ») 



(2) 



£ = (Xo/T) 1 / 2 . The first contribution corresponds to for- 
mation of a "domain wall" in the aggregate, while the 
second one corresponds to coherent tunneling transitions 
leading to re-charging of all aggregate sites. The distri- 
bution of relaxation rates can be calculated as 



As shown in 15] , the integral is dominated by the values 
N = N c , p = p c , where 

N c (r) = [ln^oT)/£] 6/5 , Pc (r) = ^(t)] 1 ^ . (4) 

The quantities iV c and p c characterize most important 
clusters among those switching during the time ~ v~ l . 
Since the number of electrons in cluster N c depends log- 
arithmically on all relevant parameters and it cannot be 
very large because other relevant parameters are expo- 
nentially sensitive to it we set N c ~ 1 following Ref. 
Substituting {4| in Eq. |3]) we obtain 



V(r) ~ ^ 



1 



T (V T) a 



(•5) 



Here P = l/(T a 3 £ 3 ), a(r) ~ r 6/5 ln 1/5 (^ r) < 1. 

Note that at 8V g 3> £h the typical distance R w 
e 2 /(n5V g ) corresponding to the polaron shift ~ 6V g pro- 
duced by a nearest neighbor turns out to be less than the 
hopping length r/j = a£, as well as a typical size ~ of 
the aggregate. Consequently one can treat the interac- 
tion between the site belonging to the percolation cluster 
and a fluctuator as a contact one. Thus, U(R) ~ e 2 / kR 
where R is the distance between the hopping site and its 
nearest neighbor belonging to the fluctuator. As result, 
the contribution of the clusters with relaxation time r to 
the distribution of polaron shifts is 



%hR 2 V{t) _ 87re 6 
d\nU/dR ~ k 3 U 3 



V(r). 



(0) 



Here we have taken into account that only aggregates 
with E < U form the polaron as well as the fact that 
each of the 2N cites of the aggregate can be coupled to 
the hopping site. The proper distribution J-(U, t) is de- 
termined by the manipulation protocol. For example, if 
the system is brought to some state at time t$ by a sud- 
den change of the gate voltage then the polaron clouds 
are formed by all the fluctuators which have changed 
their states by the observation time, t. Consequently, 

T{U,t) = Sl? r {U)dT. 

Discussion - By now we were discussing the ES 
regime of VRH. One can expect that the number of the 
metastable aggregates strongly decreases within the Mott 
regime. Indeed, aggregates are constructed from the sites 
where the intersite Coulomb energies are of the order of 
single-particle energies. This is not the case for typical 
hopping sites in the Mott VRH regime where the spread 
in the energies of the localized states is rather due to 
extrinsic disorder than to Coulomb interaction. Because 
of this spread it is less probable to find a set of sites 
forming a two-state aggregate. The addition factor en- 
tering the probability A for adding a pair of sites to an 
aggregate can be estimated as the ratio of the width of 
the Coulomb gap, Ac ~ T ' T M , to the typical hop- 



V(t) = / dNdpP(N,p)5[r(N,p)-r} . (3) ping band in ^ Mott regim6; £m = T ^ T 3/^ Here 
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Tm is the characteristic temperature of the Mott VRH, 



-(Tm/T) 1 / 4 . 



we have defined Ac as the temperature 



a oc e 

of crossover between regimes of ES and Mott VRH. The 
ratio Ac/ £m — (Ac/T) 3 / 4 is additional factor entering 
the probability A. Deeply in the Mott regime, T 3> Ac, 
this factor is small. Since Tm oc 1/po where go is the 
Mott density of states A decreases with decrease of g . 

To compare our prediction with experimental results 
of d 0, i, 0, EI M, EH obtained using InO films we 
assume that in the absence of Coulomb interaction their 
DOS, go, would be energy independent at the energies 
less or of the order of both room temperature, T r , and 
the shift in the chemical potential, S/i, due to variation 
in the gate voltage, SV g , within the dip. This DOS con- 
sists of localized and extended states split by the mobility 
edge, e m . The closer is \i to e m , the larger is the localiza- 
tion length and, consequently, the hopping conductance. 
Different samples have different go and different /x with 
respect to e m . This picture is conventional for materials 
with large amount of disorder. 

The authors of Q determine the carrier concentration 
from the Hall coefficient at T — T r . Then the found con- 
centration, n T , is just the concentration of the extended 
carriers, and one can estimate DOS as go ~ n T /T T , The 
width of the dip is related with the shift in chemical po- 
tential as SV g = goSfi/C ~ 8/j,, where G is capacitance. 
Thus the width of the dip is oc go and, for a given position 
of the chemical potential, it is correlated with n r found 
in 0]. On the other hand, the samples with the same 
go, but different positions of the chemical potential, have 
different resistance - the lower /z the lower is localization 
length and the larger is resistance. From phenomeno- 
logical point of view, the samples with the same go, but 
larger resistance can be characterized as "more dirty" . 
For density of states 50 = (4 • 10 18 cm -3 )/ (300 K) - 10 32 
cm _3 erg _1 which roughly corresponds to the threshold 
of the memory effect in 0, 0] , the Coulomb gap can be 
estimated as e c = (g e 6 / k 3 ) 1/2 ~ 3 • 10~ 15 erg. At the 
same time, T = 4 K corresponds to a crossover between 
the Mott and ES VRH regimes, and at this tempera- 
ture £~ 5 and eu ~ sc- Thus, we explain the observed 
in [B|,I3] correlation between pronounced decrease of the 
memory effect and decrease of the dip width by decrease 
of the density of the metastable aggregates in the Mott 
regime. 

According to Eq. |T]) the relative change of the con- 
ductance is x e~ 2 . Thus the shape of the dip can be cast 
in an interpolation formula 

6G {e 2 /nfPoQ{t) 



G 



(CSV g /ga) 



Q(t) 



dr (v q t)~ 



(7) 



where the limits of integration are determined by the 
manipulation protocol. The temperature dependence of 
the dip magnitude is given by the product £ _9 / 5 £/7^ oc 
T . Showing the same trend as in experiment [5| it 
is still much weaker. Assuming In vqt ~ 20 we estimate 



the height of the dip as ^^(ln i/ t)~ 6/5 ~ 0.05 that is in 
agreement with experiments. 

To analyze time dependence of the dip we take into 
account that the parameters of the system (like concen- 
tration, localization length, etc.) are somewhat different 
for different gate voltage. It is clearly demonstrated by 
the fact that the G(V) curve have a systematic slop (sub- 
tracted in course of studies of dip). Let us for simplicity 
assume that the bonding parameter £ depends on the 
gate voltage, say, through the localization length. As it 
seen from experiment, G increases with increase of con- 
centration, i.e. with V g . Consequently, we can assume 
that £ decreases with increase of V g . Correspondingly, 
the parameters of the aggregates also depend on V g and 
are different for the aggregates responsible for different 
dips in the two dip experiment. 

In the well known double-dip experiments [U the typi- 
cal relaxation time is defined as following. First, the gate 
voltage in a gated sample is rapidly changed from some 
initial to some final value V g \ (we have denoted the time 
of this variation as T m ; n ). Then it is kept constant un- 
til some time, which we will denote as i max . During this 
time the conductivity slowly (apparently logarithmically) 
decreases to some value, Go — <$G(i max , r m - m ). Then the 
gate voltage is swept to some other value, V g 2 and kept 
constant, the conductivity decreasing with time forming 
a new dip, Go — 5G(t, r m i n ). Here with a logarithmic ac- 
curacy we ascribed the same estimate r m i n for the fast 
process of switching from V g i to V g %. Let us assume that 
the shift of fj, due to the variation of V g is less than the 
typical single-particle energies of the sites forming the 
aggregates. Then at the new value of V g the aggregates 
responsible for the polaron gap at V g = V g \ stay at the 
same configuration of the occupation numbers. However 
the occupation of the sites forming the percolation cluster 
at V g — Vg\ at V g = V g 2 is completely changed provided 
that the shift of the chemical potential is larger than £/,. 
Thus the aggregates responsible for the first dip start 
to relax. However at the times t < r max the slow ag- 
gregates still preserve the configuration corresponding to 
V g = V g i; thus the first dip is partly restored if the gate 
voltage is returned to the value V g \. The depth of the 
restored dip at the time t is expected to be SG(r max ,t). 
The relaxation time, f , is defined according to equality of 
the depths of the "old" and "new" dips. This condition 
corresponds to the equality f is calculated according to 
the following procedure: 



V(t) dr 



V{t) dr . 



(8) 



Since V(r) w const x r 1 {vo T ) ° Eq. (0 then up to the 
lowest approximation in a, tq 



xAVnaxTmin- Next itera- 



tion provides the correction fi/fo — —j\n (T max /T m i n ), 
which leads to decrease of f with temperature increase. 
However, it is more sensitive to the possible dependence 
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of the parameter a on the gate voltage. The first (ini- 
tial) dip and the second one correspond to different 
gate voltages one can expect that they are formed by 
the states with different localization lengths. Conse- 
quently, the values of a are different. Denoting them 
as a\ t 2 for the first and the second dips, respectively, 
and assuming that \a% — 0:2 1 <C a we arrive at the sec- 
ond temperature-dependent correction f^/fo ~ — 
otiMAa) ln(i/ n T n )(r, Tlax /r nliT1 ) a / 2 . As follows from experi- 



ments |7J, Il7j the conductance for the second dip is larger 



than for the first one, which indicates the smaller value of 
£ for the second dip. Since a cx £~ 6 / 5 one concludes that 
a.\ > «2- Thus the correction fi is 'positive and increases 
with temperature both due to increase of a and vq. This 
trend can qualitatively explain observed weak increase of 
f with temperature. Similar conclusion can be made for 
another protocol [l8| for determining a typical relax- 
ation time. As we have seen, increase of resistance, or £ 
is correlated with decrease of the exponent a and sub- 
sequent slowing down the time evolution. This can be 
a qualitative explanation of the observed in Q slowing 
down the time evolution with increase of disorder. 

Now let us discuss an effect of external magnetic 
field |5j. One can imagine two possible mechanisms: (i) 
shrinkage of the wave function manifesting itself as a posi- 
tive addition oc H 2 to the hopping exponent [l6| ; (ii) spin 
effect related to a presence of doubly occupied centers. 
In the latter case the magnetic field align spins of singly 
occupied sites which blocks the spin-conserving hops be- 
tween the singly occupied sites. One can expect that 
at fJ,gH << Eh the latter mechanism leads to a positive 
magnetoresistance, its magnitude being proportional to 
relative contribution of the doubly occupied sites. Both 
mechanisms can be accounted for by a field-dependent 
increase of the tunneling length p entering the second 
item in the r.h.s. of Eq. (|2|). That would, in turn, lead to 
decrease of the exponent a. 

In the experiments Q , the observed magnetoresistance 
is only weakly dependent on the magnitude of resistance 
and decreases with resistance increase. This behavior 
seems to be contradictory to the wave shrinkage mecha- 
nism since in that case the magnetoresistance would dra- 
matically increase with increase of the hopping exponent 
£. In addition, the shrinkage effect is expected to be small 
for materials with small localization length. Thus, it is 
the spin mechanism that probably dominates. 

Interestingly, the double-dip memory effects are not 
observed in standard semiconductor materials. We be- 
lieve that the reason is that the ES regime of VRN in 
such materials either occurs at very low temperatures 
(less than few Kelvins) which implies weak heat with- 
drawal, or corresponds to very large resistances. Both 
of these factors seems to be disadvantageous for typical 
memory experiments. 

To conclude, our model qualitatively explains the fol- 
lowing experimentally observed features of the memory 



effect: (i) double-dip behavior of the conductance as a 
function of gate voltage; (ii) suppression of the above 
phenomenon at small carrier concentrations due to pos- 
sible crossover to the Mott regime of VRH; (iii) rather 
counter-intuitive slowing down of the time evolution (ex- 
pressed through the effective relaxation time f ) with tem- 
perature increase; (iv) slow power-law relaxation tending 
to logarithmic with increase of resistance; (v) qualitative 
dependences of memory dips on temperature and electron 
concentration; (vi) slowing down the relaxation with in- 
creasing of external magnetic field and degree of disorder 
characterized by increase of resistance. 
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